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• Digital signals, sampling, and spectra 

• Filtering 
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Spectral decomposition 



• All signalscan be 
decomposed to 
sinusoidals 






FourierTia nsform 


oo 

H(f) = \h(t)e i2xf, dt 
h(t)= j H(f)e- i2% "df 


Fouriertransorm 
Inverse Fourier transform 


Kt)= | A(/)cos(27i/ 1 + 0(f)) df Real pa rt of the signal 

-i |A(/) sm(2nft + 0(f))df imaginary part 
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Phase Magnitude 


One sided spectrum 



\ 

^ ^ l\ An 

r\ IV\^ a 

^ %A/ V 



-100 o 100 

Frequency 


For teal valued signals 

A(f) =A(-f) 

e(f) =-e(-f) 


The negative side 
does not carry 
any new informa¬ 
tion 
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Fowerspectral density 


Power spectra I density 

• Two sided: PSD = \H(f )\ 2 

• Onesided: PSD = \H(f )\ 2 + \H(-f )\ 2 

rj-i 2 

• Unite.g.- 

y Hz 

RMS spectral density 

• SD = + \H(-f )\ 2 

' Unite - g -vfc 


MJK 8/2000 



Sampling continuous signals 



Aa = amplitude resolution 
At = time resolution 


•Amplitude resolution 
•Dynamic range 
•Time resolution 
• Length in time 
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Nyquist frequent y 



Nyquist frequency 

Maximum allowed 
frequency f n =l/2f s 




Aliasing 

If f>f n it appears as 
a lowerfnequency. 
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Aliasing / Folding 



Due to the way how the 
too highfrequenciesare 
aliased, this is a Iso called 
"folding". 
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Discrete Fourierlrensfbimation 
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Leakage 



01 2345678 
Frequency [bins] 


DFTassumesa cyclic ex¬ 
tension of the data. 


Frequencies between 
bin centers "leak" also 
beyond nearest bins 
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Amplitude 


f Data Windowing 




Square 

window 



Hanning 

window 


Time 


Windowing in time domain 


Leakage 
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Frequency [bins] 


Windowing 

• Decreasesside lobes 

• Widensthe ma in lobe 
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Detrending 


Original Detrended 



• Fit a polynomial of low degree (typically a line) 
• Substractthe polynomial 
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Averaging Spectra 


PSD estimate at each bin 

•Expection value E-{P }=Mean powerin "window" 

• Standard deviation a=E-{P} 

Howto decrease the variance? 

•Increasing DFTIength does not help! 

•Shortertransforms& averaging 

• Smoothing the spectrum (combining binstogether) 


Bottom line: ^ 


Number of bins 
Number of samples 
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fr Averaging Spectra 

Hanning 



Time 


50 % overlap 
DFT lenght 


The most common way: 

• Split data span 
into overlapping 
segments 

• Apply window to 
each segment 

• Apply DFTto 
each segment 

• Square 

• Average spectra 
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FFT 

DFT 

• (Mathematical) transformation 

• Can be calculated to any data length N 

FFT 

•Algorithm to calculate DFT efficiently 
•Requireslength N =2 n 

•Numberof operations needed ~N log 2 (N) 
compared to N 2 
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Autoregressive method 


Alias 

• Maximum entropy method • All poles method 

Steps: 

• Estimate autocorrelation (non-cyclical) 

• Fit an autoregressive model to this. 

• Convert model to frequncy domain. 

Features 

• No cyclical extension 

•Can represent infinite powerdensity 

• Parametric, correct order? 



Linear fitters 


x(t) 

ax(t) 

X 1 (t)+X 2 (t) 

x(t+ t) 


N\N 





>► 


y(t) 


a y(t) 1 

> Linear 

y 1 (t)+y 2 (V J 


y(t+i) Time invariant 


vw 


Changes amplitude and 
phase but NOT frequency. 
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Frequency response 



Frequency response of a 4th o rder Butte rworth filter 





Impulse Response 
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Impulse response 


C a usa M Non-c a usa I filte is 
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